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This paper is devoted to the thermally activated dynam- 
ics of the capillary condensation. We present a simple model 
which enables us to identify the critical nucleus involved in 
the transition mechanism. This simple model is then applied 
to calculate the nucleation barrier from which we can ob- 
tain informations on the nucleation time. These results are 
compared to the numerical simulation of a Landau-Ginzburg 
model for the liquid-vapor interface combined with a Langevin 
dynamics. 
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Humidity is known to strongly affect the mechanical 
properties of many substances, like granular media or 
porous materials Jl|. Water vapor may for example con- 
dense in the pores of the medium to form liquid bridges. 
This phenomenon is the well known capillary conden- 
sation (see e.g. 0). The Laplace pressure inside such 
liquid bridges may reach a few atmospheres, and thus 
results in high adhesion forces inside the material. More 
fundamentally, capillary condensation is a confinement 
induced gas-liquid phase transition. While in the bulk 
the gas phase has a lower free energy, and is thus the 
stable phase, a liquid phase can condense in a pore when 
the liquid partially wets the solid substrate (more specif- 
ically when jsl < 7s v, where 7 is a surface tension, 
SL and SV denote the solid-liquid and the solid-vapor 
interfaces, respectively). A basic model of confinement 
is provided by the slab geometry, for which the fluid is 
confined between two parallel planar solid walls. Macro- 
scopic considerations based on this model predict a con- 
densation of the liquid phase below a critical distance H c 
between the solid surfaces satisfying the condition: 

Ap Afi ~ 2{ lsv ~ 1sl)/H c (1) 

where Ap = pl — pv is the difference between the bulk 
densities of the liquid and the gas phase, and Ap, = p sa t — 
p is the (positive) undersaturation in chemical potential 
(psat is the chemical potential at bulk coexistence) 0. 
The static part of this transition is now well understood 
and documented, both from the experimental j^Q] and 
theoretical point of view ]^|| . 

On the other hand, the dynamics of the transition have 
received very little attention. Only indirect informations 
are experimentally available in the litterature. Exper- 
imental studies of the capillary condensation using the 
Surface Force Apparatus (SFA) report strong metasta- 



bility of the gas phase when H < H c , which persists 
over macroscopic times . Furthermore the dynamics 
have been probed indirectly in experiments measuring 
the time dependent building up of cohesion inside di- 
vided materials. The latter is indeed a direct measure of 
the dynamical construction of liquid bridges between the 
grains of the medium (see ref. Q for details). From the 
theoretical side, macroscopic arguments for the topologi- 
cally equivalent drying transition predict extremely large 
time-scales for evaporation to occur 0, while the latter 
is found to occur must faster in corresponding lattice-gas 
simulations ||. 

A detailed theory for the dynamics of the capillary 
condensation providing an estimate of the condensation 
time is thus still lacking. The dynamics of the bulk gas- 
liquid transition have received on the contrary much more 
interest Away from the spinodal line, in oversatu- 
rated situations (when p > p sa t), a nucleation barrier 
can be constructed in the bulk situation, and a critical 
nucleus identified. The latter takes the form of a spherical 
droplet, with a radius Rq ~ f ApAp maximizing the 
free energy of the droplet. Since capillary condensation 
is also a first order phase transition, it should be possi- 
ble to identify a critical nucleus away from the spinodal 
line pf]| . However, the situation is more complicated in 
a confined geometry since the size of the critical nucleus 
Rq competes with other length scales, like the separation 
H between the walls. In the following, we show how to 
construct the critical nucleus for capillary condensation. 
A simplified model keeping only the main ingredients for 
capillary condensation will be first considered. The latter 
has both advantages to allow tractable calculations and 
to capture the essential features of the involved physics. 
In our case, chemical potential, total volume and temper- 
ature are fixed. Our aim is thus to find out the saddle- 
point of the grand potential of the system corresponding 
to the critical nucleus. 

To simplify further the discussion we first consider a 
two dimensionnal system and a perfect wetting situation: 
Isv — Isl + Ilv- Both assumptions shall be relaxed at 
the end. In our simplified description, we shall assume 
that the phase equilibrium is determined by macroscopic 
considerations, so that the H dependences of the surface 
tensions are totally neglected. Finally, we assume that 
the system exhibits the mirror symmetry with respect to 
the H/2 plane. Let us consider the situation in which 
planar liquid films of varying thickness e (e < H/2) de- 
velop on both solid surfaces. Following Evans et al. M, 
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the grand potential of the system may be written 



fi = - Pv Vv - p L V L + 2 1SL A + 2j LV A (2) 



where Vy (resp. Vl) is the volume of the gas (resp. 
liquid) and A is the surface area. Using Vl = 2Ae, 
Vv = A(H — 2e) and pv — Pl — ApAp,, one gets 

Auj(e) = — (fi - fi(e = 0)) = ApAp 2e (3) 

Note that in the complete wetting situation fi(e = 0) can 
be identified with fiy, the grand-potential of the system 
filled with the gas phase only. The situation e = H/2 cor- 
responds to the opposite case where the two liquid films 
merge to fill the pore. The grand potential thus exhibits a 
discontinuity at e = H/2 corresponding to the disappear- 
ance of the two liquid- vapor interfaces, and its value is 
reduced by 2jlvA. When e = H/2, expression (||) must 
then be replaced by Auj(e = H/2) = -ApAp(H c - H), 
where H c is the critical distance defined in eq. (|l|). One 
may note that the minimum of the grand potential cor- 
responds to a complete filling of the pore by the liquid 
phase when H < Hc, as expected. Up to now, liquid- 
vapor interfaces were assumed to be planar. If we allow 
deformation of the interfaces, i.e. the thickness e is now 
a function of the lateral coordinates, the corresponding 
cost has to be added to the grand potential. Due to the 
mirror symmetry assumption, the two films are identical, 
so that one finds eventually 

Afl tot = J dx { 7 Lv|Ve| 2 + A^(e)} (4) 

with AQ t ot = fi({e}) — fiy and a small slope hypothesis 
has been made. Extremalization of the grand potential 
in two dimensions leads to the following Euler-Lagrange 
equation for e(x), where x denote the lateral coordinate: 

d 2 e dAuj(e) , , 

^Ilv-t^ - — j^ 1 = (5) 
dx de 

This last equation is formally equivalent to the mechani- 
cal motion of a particle of mass 2jLVi with position e in 
the external potential — Au>(e); x plays the role of time. 
We look for solutions satisfying e = and de/dx = at 
infinity. Starting from e = 0, the "particle" is uniformly 
accelerated until it reaches e = H/2. We can choose this 
last point to fix the origin i.e. e[x = 0) = H/2. At this 
point, the discontinuity in the potential induces a spec- 
ular reflexion, similar to a collision of the particle with a 
hard wall : de/dx is therefore discontinuous and antisy- 
metric at x — 0. After a straightforward calculation the 
complete solution, depicted in fig [y, can now be obtained 
and is found to have a spatial extension x c = \/HR c , with 
R c = H c /2. Explicitly one gets e(x) — (\x\ — x c ) 2 /2R c 
for x € [—x c ; x c ] and zero otherwise. 
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FIG. 1. Picture of the critical nucleus for capillary conden- 
sation in two dimensions and perfect wetting case (6 = 0). 
The radius of curvature of the meniscus is equal to R c = H c /2, 
which is only approximatively verified within the small slope 
assumption. See text for details. 

Let us note that the cusp in the solution at x = stems 
from the discontinuity of Alu at e = H/2 resulting from 
the assumption of an infinitesimely narrow liquid-vapor 
interface. Condensation thus occurs through the exita- 
tion of short wavelength fluctuations, in agreement with 
the simulations results for the drying transition The 
corresponding energy of the nucleus (per unit length in 
the perpendicular direction) can be calculated by inte- 
grating eq. (Q), to obtain 

Afit - ^{ApA PlLV ) 1,2 H^ 2 (6) 

This energy can be identified with the energy barrier to 
overcome in order to condense the liquid phase from the 
metastable gas phase. When the energy barrier is not 
too small (compared to fegT), the time needed to con- 
dense can be estimated to be t = t cxp(Af2t//c B T), with 
To a "microscopic" time. It is easy to check that A fit 
corresponds to a saddle-point of the grand-potential. It 
is greater than both free energies of the gas and liquid 
phases. Moreover A fit is smaller than the free energy 
of any other configuration maximizing the grand poten- 
tial since it is the only solution of finite extension. Wc 
postpone the physical interpretation of the results to the 
end of the letter. We just point out that the parabolic 
solution obtained above is the small slope approximation 
to the circle with radius of curvature R c . 

In order to verify the previous results, we have con- 
ducted numerical simulations of the capillary conden- 
sation in two dimensions. We start with a Landau- 
Ginzburg model for the grand potential of the system 
confined between two wall. In terms of the local density 
p(r), we write the "excess" part of the grand potential 
fl ex — fl+P sat V, where P sa t is the pressure of the system 
at coexistence, as 

fi e * = J dv{^p\ 2 + W{p) + {Ap + V ext {z))p) (7) 
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In this equation, m is a phenomenological parameter; 
Vext(z) is the confining external potential, which we took 
for each wall as V ext {z) = —e(a/(Az + er)) 3 , with Az 
the distance to the corresponding wall; e and a have 
the dimensions of an energy and a distance. W(p) can 
be interpreted as the negative of the "excess" pressure 
Msatp- f(p)-Psat, with f(p) the free-energy density 0. 
As usually done, we assume a phenomenological double 
well form for W(p) : W(p) = a(p — pv) 2 {p — Pl) 2 , where 
a is a phenomenological parameter jl2]. The system is 
then driven by a non-conserved Langevin equation for p : 

where T is a phenomenological friction coefficient and 
77 is a Gaussian noise field related to T through the 
fluctuation-dissipation relationship Jl3[ |. An equivalent 
model has been successfully used for the (bulk) classic 
nucleation problem [jlij. Physically, the non-conserved 
dynamics assume an infinitely fast transport of matter 
in the system, which is justified in view of the time-scale 
involved for condensation. We solved (||) by numerical 
integration using standard methods, identical to those 
of ref. The units of energy, length are such that 

(7 = 6 = 1 . Time is in unit of to = (recT 2 )^ 1 with T = -|. 
In these units, we took m = 1.66, a = 3.33, pl = 1, 
pv — 0.1. Typical values of the chemical potential and 
temperature are A/i ~ 0.016, T ~ 0.06 (which is roughly 
half the energy barrier between vapor and liquid with the 
form for W(p) used in our model). Periodic boundary 
condition with periodicity L x were applied in the lateral 
direction. 



that the amplitude of capillary waves increases with the 
lateral dimension of the system for free interfaces |Q| . In 
our case however, the long-range effects of the fluctua- 
tions of the liquid film are expected to be screened due 
to the presence of the external potential. Moreover, as 
predicted by the model, nucleation should occur via the 
excitation of localized fluctuations. The observed insen- 
sitivity of the results with respect to finite size effects is 
then an encouraging feature for the model presented in 
this letter. 

The simulated system is initially a gas state filling the 
whole pore, and its evolution is described by eq. (||). A 
typical evolution of the mean density in the slit p(t) is 
plotted on fig |[ An average over different realizations 
(from 10 to 30) is next performed to get an averaged 
time-dependent density p(t). As expected due to the 
long range nature of the external potential a thick liquid 
film of thickness I rapidly forms on both walls on a short 
time scale ri (I ~ 3.8tr and n « 5io in our case). The 
first step of the dynamics is thus the wetting of the solid 
substrates. This process is not thermally activated. In a 
second stage, fluctuations of the interfaces around their 
mean value I induce after a while a sudden coalescence of 
the films (see fig. ||). This second process has a charac- 
teristic time r. It is numerically convenient to define the 
total coalescence time, n + t, as the time for the average 
density in the slab between the two wetting films to reach 
(pv + Pl)/2 H|, which corresponds in our case to the 
condition p(r + t\) « 0.8. The physical results do not 
depend anyway on the precise definition of t. In fig. |^, 
we plot the variation of ln(r) as a function of the inverse 
temperature 1/T. 
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FIG. 2. Averaged density as a function of time t (in units 
of to) for a few realizations of the noise (H = 13, Ap = 0.016 
and T = 0.07). The dashed line is the average over all the 
realizations, p(t). 

Typically L x ~ 2H was used, but we have checked that 
increasing L x up to 20H does not affect the results for 
the activation dynamics. We emphasize that this lack 
of sensitivity is not an obvious result since it is known 
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FIG. 3. Logarithm of condensation time as a function of 
the inverse temperature (A/i = 0.016, H = 13). The dashed 
line is least-square fit of the datas. 

As expected, far from the spinodal (i.e., for large 
enough H, iJ>,3£), r is found to obey an Arrhenius law 
t = tq exp(Afl' ! /ksT), where AW is identified as the 
energy barrier for nucleation. We now focus on the vari- 
ations of ln(r) as a function of L x and H, which one 
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assumes to be mainly controlled by the variation of Afit. 
First, as already noticed above, we found no variation 
of ln(r) as a function of L x , in agreement with our pre- 
diction of a localized critical nucleus. The dependence 
on H (Afj, being fixed) is plotted on fig. [|. The pre- 
vious model predicts a H 3 / 2 dependence. However the 
long range of the external potential (of the van der Waals 
type), produces thick wetting films whose thickness has 
to be substracted from H . A more careful analysis of 
the critical nucleus shows in fact that the total effective 
thickness of the films has to be replaced by M (instead 
of 2£), in order to take correctly into account the long 
range character of the external potential. This result is in 
agreement with other theoretical and experimental find- 
ings for capillary condensation in the presence of van der 
Waals forces f|Jlj|. As seen on fig. |[ a good agreement 
with the theoretical prediction is found. Dependence on 
the other parameters (A/i, e, ...) shall be discussed in a 
longer version of this paper. 



Af2 f =7 L yH c {a-sinacos(a;-t-260} (9) 

with a defined through cos(a + 9) = cos(#) — H/H c . 
The previous result, eq. (||), is recovered in the lim- 
iting case 9 — 0, H -C H c . In three dimensions, the 
nucleus takes the form of a liquid bridge of finite lat- 
eral extension connecting the two solid surfaces, due in 
particular to the supplementary (negative) axisymmetric 
curvature. When H is small compared to H Cl this pre- 
dicts Afit (x -fLvH 2 + C(A/i), but the full dependence 
on Afj, for any H needs a numerical resolution. This will 
be done in a forthcoming paper together with the corre- 
sponding simulations. 
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ENS-Lyon, and the MENRT under contract 98B0316. 
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FIG. 4. Logarithm of condensation time as a function 
of the "effective" width of the slab H — 31 for fixed 
Afi — 0.016. The dashed line is the theoretical prediction 
ln(r) = ln(ro) + ct(H — 3£) 3 ^ 2 . The two parameters ln(ro) and 
a have been obtained from a least-square fit of the datas in a 
ln(r) versus (H - 3f) 3/2 plot. 

The simulation results in the 2D perfect wetting case 
are thus in agreement with the theoretical predictions. 
Now it is possible to relax the assumptions done in the 
presentation above, i.e. perfect wetting, two dimension- 
nal system, small slopes. More generally, one may realize 
that the extremalization of the grand potential leads in 
one hand to the usual Laplace equation, relating the local 
curvature k to the pressure drop jlv k = Ap ~ A/iAp; 
and in the other hand, it fixes the contact angle of the 
meniscus on the solid substrate according to Young's law 
7l\/cos# = j sv — jsl- In our case however, the corre- 
sponding nucleus corresponds to a maximum of the grand 
potential. In two dimension the general solution is of the 
same geometrical form as the one obtained above and 
pictured in fig |l|, with a corresponding nucleation energy 
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